LOWER BOUNDS FOR POLYNOMIALS USING GEOMETRIC 

PROGRAMMING 



MEHDI GHASEMI AND MURRAY MARSHALL 

Abstract. We make use of a result of Hurwitz and Reznick [8] [19] . and a 
consequence of this result due to Fidalgo and Kovacec [5] , to determine a new 
sufficient condition for a polynomial / £ M[Xi, . . . , X n ] of even degree to be a 
sum of squares. This result generalizes a result of Lasserre in [10] and a result 
of Fidalgo and Kovacec in [5] , and it also generalizes the improvements of these 
results given in [6]. We apply this result to obtain a new lower bound f gp for 
/, and we explain how f gp can be computed using geometric programming. 
The lower bound f gp is generally not as good as the lower bound f sos intro- 
duced by Lasserre |11| and Parrilo and Sturmfels |15| , which is computed using 
semidefinite programming, but a run time comparison shows that, in practice, 
the computation of f gp is much faster. The computation is simplest when the 
highest degree term of / has the form JZILi "l^'f : a i > 0, i = 1, . . . , n. The 
lower bounds for / established in [6] are obtained by evaluating the objective 
function of the geometric program at the appropriate feasible points. 



1. Introduction 

Fix a non-constant polynomial / G M[X] = R[Xi, ■ ■ ■ ,X n ], where n > 1 is an 
integer number, and let /* be the global minimum of /, defined by 

/* := inf{/(o) : a € R n }. 

We say / is positive semidefinite (PSD) if f(a) > Va G R". Clearly 

inf{/(a) : a G R"} = sup{r G K : / - r is PSD}, 

so finding /* reduces to determining when / — r is PSD. 

Suppose that deg(/) = m and decompose / as / = /o H + f m where /, is a 

form with deg(/i) = i, i = 0, . . . , to. This decomposition is called the homogeneous 
decomposition of /. A necessary condition for /* ^ -co is that f m is PSD (hence 
to is even). A form g G M.[X] is said to be positive definite (PD) if g(a) > for all 
a G R™, a ^ 0. A sufficient condition for /» ^ — oo is that f m is PD [14]. 

It is known that deciding when a polynomial is PSD is NP-hard [TJ Theorem 1.1]. 
Deciding when a polynomial is a sums of squares (SOS) is much easier. Actually, 
there is a polynomial time method, known as semidefinite programming (SDP), 
which can be used to decide when a polynomial / G K[A] is SOS [TT] [TH). Note 
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that any SOS polynomial is obviously PSD, so it is natural to ask if the converse is 
true, i.e. is every PSD polynomial SOS? This question first appeared in Minkowski's 
thesis and he guessed that in general the answer is NO. Later, in [TJ, Hilbert gave 
a complete answer to this question, see [2j Section 6.6]. Let us denote the cone of 
PSD forms of degree 2d in n variables by P2d,n and the cone of SOS forms of degree 
2d in n variables by E2d,n- Hilbert proved that Pid,n = ^id,n if and only if (n < 2) 
or (d = 1) or (n = 3 and d = 2). 

Let J2 MX? denote the cone of all SOS polynomials in R\X] and, for / G R\X], 



Since SOS implies PSD, f sos < /*. Moreover, if f sos ^ — oo then f sos can be 
computed in polynomial time, as close as desired, using SDP pj] [15]. We denote 
by P 2 ° d n and T^ d n , the interior of P2d,n and £2d,n in the vector space of forms of 
degree 2d in R[A], equipped with the euclidean topology. A necessary condition 
for f sos — oo is that fid G ^2d,n- A sufficient condition for f sos ^ — oo is that 
hd G „ [El Proposition. 5.1]. 

In Section [3 we recall the Hurwitz-Reznick result (Theorem l2.ip and a corollary 
of the Hurwitz-Reznick result due to Fidalgo and Kovacek (Corollary 12. 2|) . For the 
convenience of the reader we include proofs of these results. Using the latter result, 
we determine a sufficient condition, in terms of the coefficients, for a form / of 
degree 2d to be SOS (Theorem I2.3[) . We explain how Theorem 12.31 can be applied 
to derive various concrete criteria, in terms of the coefficients, for a form to be SOS, 
including results proved earlier by Lasserre pjQl Theorem 3], Fidalgo and Kovacec 
Theorem 4.3], and Ghasemi and Marshall j6j Section 2]. 

In Section [3l we use Theorem 12.31 to establish a new lower bound f gp for / and 
we explain how f gp can be computed using geometric programming. An advantage 
of the method is that solving a geometric program is almost as fast as solving a 
linear program. Although the lower bound found by this method is typically not 
as good as the lower bound found using SDP, a practical comparison shows that 
the computation is much faster, and larger problems can be handled. 

In Scction|4]we explain how results in Section[3]imply and improve on the results 
in [6j Section 3]. 

In this paper we denote by N the set of nonnegative integers {0, 1,2,.. .}. For 
A = (Xi, . . . ,X n ), a = (ai,...,a„) e M™ and a = {oi\, . . . ,a n ) G N", define 
X_ a := A" 1 • • • A"", |a| :=£*! + ••■ + a n and a a :— a" 1 • • • a" n with the convention 
0° = 1. Clearly, using these notations, every polynomial / £ M[A] can be written 
as /(A) = X)q6N" /a^": where f a G M and f a — 0, except for finitely many 
a. Assume now that / is non-constant and has even degree. Let f2(/) = {a G 
N" : f a + 0} \ {0, 2de x , . . . , 2de n }, where 2d = deg(/), = (6 lU . ..,S in ), and 



define 



fsos '■= sup{r G R : /-r£^K[I] 2 }. 
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We denote /o and /2dei by /o and fid,i for short. Thus / has the form 

n 

(i) f = fo+ E /«* q +Em^ 

oen(/) t=i 

Let A(/) = {a e Sl(f) : / Q X" is not a square in = {a e : either f a < 

or a, is odd for some 1 < i < n}. Since the polynomial / is usually fixed, we will 
often denote and A(/) just by f2 and A for short. 

Let f(X, Y) = Y 2d f(^, . . . , From © it is clear that 

QGf2 i=l 

is a form of degree 2d, called the homogenization of /. We have the following 
well-known result: 

Proposition 1.1. / is PSD if and only if f is PSD. f is SOS if and only if f is 
SOS. 

Proof. See QH Proposition 1.2.4]. □ 



2. Sufficient conditions for a form to be SOS 
We recall the following result, due to Hurwitz and Rcznick. 

Theorem 2.1 (Hurwitz-Reznick). Suppose p(X) = J27=i ^ x f d ~ ZdX? 1 ■ ■ ■ X%™ , 
where a — (a±, . . . , a n ) 6 N™ ; \a\ = 2d. Then p is SOBS. 

Here, SOBS is shorthand for a sum of binomial squares, i.e., a sum of squares of 
the form {aX a ~ bXf) 2 

In his 1891 paper [8], Hurwitz uses symmetric polynomials in X\, . . . , X2d to give 
an explicit representation of X^iLi — 2d\^ i=1 Xi as a sum of squares. Theorem 
12.11 can be deduced from this representation. Theorem 12.11 can also be deduced 
from results in [TBI [19], specially, from (T9l Theorems 2.2 and 4.4]. Here is another 
proof. 

Proof. By induction on n. If n = 1 then p = and the result is clear. Assume now 
that n > 2. We can assume each cti is strictly positive, otherwise, we reduce to a 
case with at most n — 1 variables. 

Case 1: Suppose that there exist 1 < ii,ia < n, such that %\ ^ 12, with < d 
and cti 2 < d. Decompose a — (a±, . . . , a n ) as a = f3 + 7 where /3, 7 G N™, = 0, 
7i2 = and \/3\ = \j\ = d. Then 



(xf - 2C) 2 = x 20 - 2x! i xj + x 21 = x_ 213 - 2X a + x 21 , 
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therefore, 

n 

p(2Q = aiX ? d ~ 2dX - a 

i=l 
n 

; 2/3iX? d - 2dX 2 ^ 

+ 5 \J2 2 ^ x ? d - 2d x 2 ^ + - 2C) 2 - 

Each term is SOBS, by induction hypothesis. 

Case 2: Suppose we are not in Case 1. Since there is at most one i satisfying 
eti > d, it follows that n = 2, so p(X) = ai X 2d + a 2 X% d - 2dX^X%\ We know 
that p > on R 2 , by the arithmetic-geometric inequality. Since n = 2 and p is 
homogeneous, it follows that p is SOS. 

Showing p is SOBS, requires more work. Denote by AGl(2,d) the set of all 
homogeneous polynomials of the formp = ctiX 2d + a2X 2d — 2dX" 1 X^ 2 , ol\,oli £N 
and ot\ + ai — 2d. This set is finite. If ct\ = or oti = 2d then p = which is 
trivially SOBS. If a x = a 2 = d then pQQ = d(Xf - X 2 d ) 2 , which is also SOBS. 
Suppose now that < u\ < 2d, a% ^ d and a± > a 2 (The argument for o.\ < a 2 is 
similar). Decompose a — (ai,a 2 ) as a = j3 + 7, /3 = (d, 0) and 7 = («i — d, a 2 ). 
Expand p as in the proof of Case 1 to obtain 

p(X) = 2 h X i d ~ 2dX 2 ^ 2 ^ X i d ~ 2dX 2 ^ + d(X? - X~<) 2 . 

Observe that £? =1 WiX 2d - 2dX 2p = 0. 

Thus p = |pi + d(X^ - X 7 ) 2 , where p x = £? =1 2 7i X? d - 2dX 27 . If Pl is 
SOBS then p is also SOBS. If pi is not SOBS then we can repeat to get p\ = 
\p2 + d(Xf — X 1 ) 2 . Continuing in this way we get a sequence p = po,pi,p 2 , • • ■ 
with each pi an element of the finite set AGI(2, d), so Pi = Pj for some i < j. Since 
Pi = 2 l ~ J pj+ a sum of binomial squares, this implies pi is SOBS and hence that p 
is SOBS. □ 

In [5], Fidalgo and Kovacec prove the following result, which is a corollary of the 
Hurwitz-Reznick result. 

Corollary 2.2 (Fidalgo-Kovacek) . For a form p{X) = YJL\^ X f d - M=£° such 
that a € N™, \a\ — 2d, ft > for i = 1, • ■ • ,n, and fi > if all are even, the 
following are equivalent: 

(1) p is PSD. 

(2) ^ d m=^r<(2d) 2d uti(3r- 

(3) p is SOBS. 

(4) p is SOS. 
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Proof. See Theorem 2.3]. (3) =>■ (4) and (4) =>• (1) are trivial, so it suffices 
to show (1) => (2) and (2) =>• (3). If some a.i is odd then, making the change of 
variables 1$ = — Xj, lj = Xj for j ^ i, fi gets replaced by — /i. In this way, we can 
assume /i > 0. If some a.; is zero, set Xj = and proceed by induction on n. In 
this way, we can assume on > 0, i = 1, . . . , n. If /i = the result is trivially true, 
so we can assume fj, > 0. If some ft is zero, then (2) fails. Setting Xj — 1 for j ^ i, 
and letting Xj — > oo, we see that (1) also fails. Thus the result is trivially true in 
this case. Thus we can assume ft > 0, i = 1, . . . , n. 

(1) (2). Assume (1), so p{x) > for all x e M™. Taking 

*:=((f) 1/M , ■•■>(f L ) 1/2d ), 

Pi Pn 

we see that 

n n n 

s° A t nr=i(fr) Q * /2<i - 2d - This p r ° ves (2)- 

(2) => (3). Make a change of variables X l = (f^) 1/2d Yi, i = l,...n. Let 
Mi := vT\LiW a ' /2d so ' b y ( 2 ). W ^ 2d ' Le -' If ^ L Then 

n n ~j n 

pQQ = E - ^ = SE - !) + E ^ - 2 « 

i=i i=i t=i 

which is SOBS, by the Hurwitz-Reznick result. This proves (3). □ 

Next, we prove our main new result of this section, which gives a sufficient 
condition on the coefficients for a polynomial to be a sum of squares. 

Theorem 2.3. Suppose f is a form of degree 2d. A sufficient condition for f to 
be SOBS is that there exist nonnegative real numbers a a> i for a € A, i = 1, . . . , n 
such that 

(1) Vae A {2df d a a a = fl d a a . 

(2) f 2 d,i > E«ga a <*,i> i=l,..-,n. 
Here, a a := (a a> i, . . . , a Qi „). 

Proof. Suppose that such real numbers exist. Then condition (1) together with 
Corollary O implies that £" =1 a a<i Xf d + f a X a is SOBS for each a € A, so 

n 

E(E a -)^ 2d +E/^ Q 

i=l q£A a£A 

is SOBS. Combining with (2), it follows that £™ =1 f 2 daXf d + £ QeA f*2C is SOBS. 
Since each f a X_ a for a 6 D, \ A is a square, this implies /(X) is SOBS. □ 

Remark 2.4. (i) From condition (1) of Theorem l2.3l we see that a a ,i = =>• on = 0. 
(ii) Let a be an array of real numbers satisfying the conditions of Theorem l2.3[ and 
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define the array a* — (a* t ) by 



a a ,i if oti ^ 
if ou = 0. 



Then a* also satisfies the conditions of Theorem 12.31 Thus we are free to require 
the converse condition on — => a a j = too, if we want. 



We mention some corollaries of Theorem l2.3l Corollaries l2.5l and l2.6l were known 
earlier. Corollary 12.71 is an improved version of Corollary 12.61 Corollary 12.91 is a 
new result. 

Corollary 2.5. See [H Theorem 3] and H Theorem 2.2]. For any polynomial 
f e U\X] of degree 2d, if 

(LI) /„> £ l/al 2 ^ and (L2) f 2dA > £ \f a \%, i = l,...,n, 
then f is a sum of squares. 

Proof. Apply Theorem 12.31 to the homogenization f(X, Y) of /, taking a a ^ = 
m> * = •"■> • • • i n an< ^ a a,Y — \fa\ 2d ^ 0T eacn a S A. For a € A, 

(2d) 2 ^ = (2d) M ( l/ " l( ^-" ,|) ) M " H nr = i( J ^r i ) a ' 

= (2d) 2d |/ Q | 2d -l Q l(2d- |a|) 2d -l Q l|/ Q |l«la a (2d)- 2d 
= \f a \ 2d a a (2d~\a\) 2d -^. 

So, IQTl) holds. (LI) and (L2) imply E2K2), therefore, by Theorem El / and 
hence / is SOBS. □ 

Corollary 2.6. See [5, Theorem 4.3] and jU Theorem 2.3]. Suppose f E R[X] is 
a form of degree 2d and 



min / 2d ,i>i^|Ma Q )^. 
—i,...,n zd * — ; 



Proo/. Apply Theorem El with a 0)i = [fa] 9 ^-, VaeA,t = l...,n. □ 
Corollary 2.7. Suppose f is a form of degree 2d, fid.i > 0, i — 1, . . . , n and 



\f a \a a ' 2d 

Then f is SOBS 



E \Ja\ u - ' < 1 

OrlY\ n f a i/ 2d ~ 

qGA zu lli=l J2d,i 



Proof. Apply Theorem El with a a , = l/ ° l "° /2d / 2 f^ ■ □ 

2dn? =1 / 2d j j 
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Remark 2.8. Corollary 12. 71 is an improved version of Corollarv l2.6l This requires 
some explanation. Suppose that f2d,i > Sqga \fa\& a ^ 2d , i = l,..-,n. Let 
hd,io : = min {/2d,i :i = l,...,n}. Then 



2d,io i 



and 



\f a \a a/2d = l_ |/>° /2rf / 2d , 
1 \f a \a a ' 2d 



n rOLijld 



We note yet another sufficient condition for SOS-ness. 
Corollary 2.9. Let f E R\X] be a form of degree 2d. If 

(if | \ 2(1/ ain a 
■gj) , =l,-.-,n 

iften / is SOBS. Here n a := \{i : cti ^ 0}|. 
Proof. Apply Theorem 12.31 with 

/If |\ 2d /"i"c 

if a* = 0. 

□ 

The following example shows that the above corollaries are not as strong, either 
individually or collectively, as Theorem [273] itself . 

Example 2.10. Let f(X, Y, Z) — X 6 +Y 6 +Z 6 ~5X—4:Y-Z+8. Coro]laiy|27Jdoes 
not apply to /, actually (LI) fails. Also, Corollaries 12.61 [2771 and [2791 do not apply 
to /, the homogenization of /. We try to apply Theorem 12.31 Let a% = (1,0,0,5), 
Q'2 = (0, 1, 0, 5) and 013 = (0, 0, 1, 5), then A = {a.\, 012, 013}. Denote a ai j by ay, we 
have to find positive reals 011,022,033,014,024,034 such that the followings hold: 

6 6 onof 4 = 5 6 5 5 , l>Ou, 
6 6 o 22 o| 4 = 4 6 5 5 , l>o 22 , 
6 6 a 3 3a| 4 = 5 5 , 1 > a 33 , 

8 > 014 + a 2 4 + a 34 . 

Take on = 022 = 033 = 1 and solve equations on above set of conditions, we get 
Oi4 + 024 + O34 « 7.674 < 8. This implies that / and hence / is SOBS. 
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3. Application to global optimization 

Let / € M.[X] be a non-constant polynomial of degree 2d. Recall that f sos 
denotes the supremum of all real numbers r such that /-r 6 ^] R[2C| 2 , /* denotes 
the infimum of the set {/(a) : a G R n }, and f sos < /*. 

Suppose / denotes the array of coefficients of non-constant terms of / and fo 
denotes the constant term of /. Suppose $(/, fo) is a formula in terms of coefficients 
of / such that $(/, fo) implies / is SOS. For such a criterion we have 

Vr ($(£/„ -r)-> r < f aoa ), 

so /$ := sup{r € R : $(/,/o — r)} is a lower bound for f sos and, consequently, for 
/* . In this section we develop this idea, using Theorem 12.31 to find a new lower 
bound for /. 

Theorem 3.1. Let f be a non- constant polynomial of degree 2d and r £R. Suppose 
there exist nonnegative real numbers a a i , a £ A, i = 1, . . . , n, a a< i = iff oti = 0, 
such that 

(1) {2d) 2d a% = \f a \ 2d a a for each a e A such that \a\ = 2d, 

(2) fid,% > E qG a a <« for i = 1, . . . ,n, and 

[fa] 



(3) /o-r>E« 6 A<«(2d- 



(2d) 2d a 

Then f -r is SOBS. Here A <2d := {a e A : |a| < 2d}. 



Proof. Apply Theorem 12.31 to g := / — r, the homogenization of / — r. Since 
/ = /o + E?=i f2d, l Xf d + J2 a£ n f a 2C, it follows that g = (f -r)Y 2d + J^ =1 Xf d + 
J2 a en f a X a Y 2d ~\ a \. We know f-r is SOBS if and only if g is SOBS. The sufficient 
condition for g to be SOBS given by Theorem 12.31 is that there exist non-negative 
real numbers a a ^ and a Q y, a a ,i = iff a* = 0, a Q y = iff |a| = 2d such that 

(1) ' V a e A (2d) 2d a>' d v : |Q ' = |/ Q | 2d a a (2d - |a|) 2d -H, and 

(2) ' /2d,i > Eqga i = 1, • ■ • , n and fo~r> J2 ae A a *,Y- 
Solving (1)' for a a ^y yields 



a a y = [2d — \c 



{2d) 2d a« 

if | of | < 2d. Take a a y = if \a\ — 2d. Conversely, defining a a y in this way, for 
each a S A, it is easy to see that (1), (2), and (3) imply (1)' and (2)'. □ 

Definition 3.2. For a non-constant polynomial / of degree 2d we define 

f gp := sup{r £ R : 3a ai i e R-°, a 6 A, i = 1, . . . , n, a a< i = iff o.\ = 

satisfying conditions (1), (2) and (3) of Theorem 13. II} . 

It follows, as a consequence of Theorem 13. 1[ that f gp < f sos . 

Example 3.3. Let f(X, Y) = X 4 + Y 4 - X 2 Y 2 + X + Y. Here, A = {a u a 2 , a 3 }, 
where a-y — (1,0), «2 = (0,1) and a?, — (2,2). We are looking for non-negative 
reals a t j, i — 1,2, 3, j = 1,2 satisfying an + a 2 \ + a 31 < 1, a 12 + a 22 + 032 < 1, 
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031032 = \. Taking an = a 2 2 = a 3 i = a 32 = \, a%2 = &2\ = 0, we see that 
fgp > -5^3 ■ Taking X = Y = we see that /* < f{—^, -^) = -5273 ■ 

Since f gp < f sos < /», it follows that / 9P = / sos = /» = -^j. 

Remark 3.4. If |Q| = 1 then /* = / sos = / ffP . 

Proof. Say ft = {a}, so / = £™ =0 f2d,iX? d + f + f a X a . We know / gp < f sos < 
so it suffices to show that, for each real number r, /» > r =>■ / 9P > r. Fix r and 
assume /» > r. We want to show / sp > r, i.e., that r satisfies the constrains of 
Theorem 13. II Let g denote the homogenization of / — r, i.e., g — Y%=i hd,iX 2d + 
(fa - r)Y 2d + f a X a Y 2d -\ a \. Thus g is PSD. This implies, in particular, that 
f2d,i > 0, i = 1, . . . , n and /o > r. There are two cases to consider. 

Case 1. Suppose f a > and all a* are even. Then a ^ A, so A = 0. In this case 
r satisfies trivially the constraints of Theorem 13. 11 so f gp > r. 

Case 2. Suppose either / Q < or not all of the a, are even. Then a G A, i.e., 
A = Q = {a}. In this case, applying Corollary |2.2[ we deduce that 

n 

(2) f a d a a (2d-\a\) 2d -^ < (2d) 2d JJ - r) M ~H 

i=l 

There are two subcases to consider. If |a| < 2d then r satisfies the constraints of 
Theorem 13. li taking 

_ J /2d,i if o» # 
fl "' 1 -\ if Qj = 0. 

If |a| = 2d then © reduces to /^ d a a < (2d) 2d H" =1 In this case ' r satisnes 

the constraints of Theorem 13 - H taking 

_ J s/ad.i if a t 7^ 
" \ if a* = 0. 

where 

|/a| M a a 



(2rf) 2d nr= 1 / 2 % 



□ 

If /2d,i > 0, i = 1, ... ,n then computation of f gp is a geometric programming 
problem. We explain this now. 

Definition 3.5. (geometric program) 
(1) A function / : M^ -)• E of the form 

where c > 0, ^ £ K and x = (x\, . . . , x n ) is called a monomial function. A sum of 
monomial functions, i.e., a function of the form 



k 

(x) = Yc iX ? 
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where cj > for i = 1, . . . , k, is called a posynomial function. 

(2) An optimization problem of the form 

{Minimize <fio (x) 
Subject to 4>i{x) < 1, i = l,...,m and ^i(x) = 1, i = 1,.. . ,p 

where 0o, ■ • ■ , 0m are posynomials and ipi, . . . , ip p are monomial functions, is called 
a geometric program. 

See [21 Section 4.5] or [THl Section 5.3] for detail on geometric programs. 

Corollary 3.6. Let f be a non-constant polynomial of degree 2d with f2d,i > 0, 
i = 1, . . . , n. TTien / sp = fo~ m* where m* is the output of the geometric program 



Minimize J2aeA< 2d i^d _ M) 



Subject to J2aeA jffl ^ h « = 1, • • • ,n and [^| 2d °g =1, a e A, |a| = 2d. 

The variables in the program are the a a ^, a £ A, i = a^O, the under- 

standing being that a a ,i = iff a, = 0. 

Proof. f gp = /o — to* is immediate from the definition of / sp . Observe that 

i 

2d 



<Ma):= J2 ( 2d -|"D 

aGA,|a|<2<i 



and 0j(a) := X) Q eA Tfr 7 ' * = ■"•> • • ■ > n are posynomials in the variables a a ,i, and 
"00(0) := |^ d |2d°Q , oi S A, |a| = 2c? are monomial functions in the variables a Qj i. □ 

Addendum: If either f 2 d,% < for some i or f 2 d,i = and a* ^ for some i 
and some a then / gp = — oo. In all remaining cases, after deleting the columns of 
the array {a a ^) corresponding to the indices i such that f 2 d,i — 0, we are reduced 
to the case where fid,% > for all i, i.e., we can apply geometric programming to 
compute fgp. 

A special case occurs when f 2 d,i > 0, for i = 1, . . . , n and {a <E A : \a\ = 2d} = 
0. In this case, the equality constraints in the computation of m* are vacuous and 
the feasibility set is always non-empty, so f gp 7^—00. 

Corollary 3.7. If \a\ < 2d for each a G A and f 2 d,i > for i = 1, ...,n, then 
fgp 7^ — 00 and fg P = fo — to* where m* is the output of the geometric program 



Minimize S a eA( 2 ^~ l ( 

Subject to J2a£A a a,i < fcd,i, i = !,-■■ ,n. 
Proof. Immediate from Corollary 13.61 □ 
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Example 3.8. (1) Let / be the polynomial of Example 12.101 Then f gp = f sos 
i. w 0.3265. 

(2) For g(X, Y, Z) = X 6 + Y 6 + Z e + X 2 YZ 2 - X 4 - Y 4 - Z A - YZ 3 - XY 2 + 2, 
g* « 0.667, and g gp = g sos w -1.6728. 

(3) For h(X,Y,Z) = g{X,Y,Z) + X 2 , we have h gp » -1.6728 < h sos w -0.5028 
and h* « 0.839. 

To compare the running time efficiency of computation of f sos using semidefinite 
programming with computation of f gp using geometric programming, we set up a 
test to keep track of the running times. All the polynomials were taken randomly of 

the form X\ d H V X 2d + g(X) where g G R[X] is of degree < 2d- 1. In each case 

the computation is done for 50 polynomial with coefficients uniformly distributed 
on a certain symmetric interval, using SosTools and GPposy for Matla^. 

Although, sometimes there is a large gap between f sos and f gp , the running time 
tables show that computation of f gp is much faster than f sos . 



Table 1. Average running time (seconds) to calculate /. 



n\2d 


4 


6 


8 


10 


12 


3 


0.73 


1 


1.66 


2.9 


6.38 


4 


0.98 


1.8 


5.7 


25.5 




5 


1.43 


4.13 


44.6 






6 


1.59 


13.24 


573 







Table 2. Average running time (seconds) to calculate /, 



n\2d 


4 


6 


8 


10 


12 


3 


0.08 


0.08 


0.12 


0.28 


0.36 


4 


0.08 


0.13 


0.3 


0.76 


2.23 


5 


0.08 


0.25 


0.8 


3.42 




6 


0.09 


0.37 


2.2 







Example 3.9. Let f(X, Y, Z) = X i0 + Y 40 + Z 40 - XYZ. According to Remark 
13.41 /* = / sos = fgp. The running time for computing f gp w —0.686 using geometric 
programming was 0.18 seconds, but when we attempted to compute / sos directly, 
using semidefinite programming, the machine ran out of memory and halted, after 
about 4 hours. 



hardware and Software specifications. Processor: Intel® Core™2 Duo CPU P8400 @ 
2.26GHz, Memory: 2 GB, OS: Ubuntu 10.04-32 bit, Matlab: 7.9.0.529 (R2009b) 
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4. Explicit lower bounds 



We explain how the lower bounds for / established in [6, Section 3] can be 
obtained by evaluating the objective function of the geometric program in Corollary 
13.71 at suitably chosen feasible points. 

Recall that for a (univariate) polynomial of the form p(t) = t™ — Y17=o a ^*> 
where each ai is nonnegative and at least one is nonzero, C(p) denotes the 
unique positive root of p [T7| Theorem 1.1.3]. See g], Ex. 4.6.2: 20] or |5J 
Proposition 1.2] for more details. 

Corollary 4.1. // |a| < 2d for each a G A and fid-i > for i = 1, ...,n, then 
f gp > r L , where 

r L := U-hT, a eAW-\a\)\f a \kW{f^ 
k > .m^ n C(t 2d -^Y, ae A^\fM2 d ft ]a] )- 

Here,/ 2 7:=nr= 1 / 2d > 
Proof. For each a £ A and i = 1, • • • , n. Let 

By definition of fc, for each i, ± J2aeA &t\fa\{f2d,i)~^ k^ < k 2d , hence 



2dk 2d -\ a \ 



This shows that the array (a a i : a € A,i = 1, •■• ,n) is a feasible point for 
the geometric program in the statement of Corollary 13.71 Plugging this into the 
objective function of the program yields 



E Q6 A(2d-M 



(it) Ila^O (^~) 
(it) Ila^O (a.l/flfciV^ (fod.i) 1 ^ *) 

(!tf iWo (^^ H (/-)^) ar 



= iiE a eA(2d~\a\)\f a \kW(f- d a )^, 
so r L - /o - j- d T. a eA^d-\a\)\f a \k^\{f- d a )h < f ; 



up- 



Corollary 4.2. // |a| < 2d for each a 6 A and /2<i,i > /or i = 1, . . 

/ OT > fFK; where r FK := f - k 2d , k > C(t 2d - E^" 1 bit 1 ), 

a£ A, |a|— * 



□ 



Proof. Define 



{2d-\a\) ^ —(a f 2d ) ' f 2d ,ik [ 1 
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Note that Y^a=i 1 ^ik* — k 2d and, for each i = 1, . . . , n, 

J2a a ,i = E( 2d -H)^W( aa/ 2"/) A /2^ 



-2,/ 



aeA 



aeA 
2d-l 



{2d -j)*&-\M (a°7 2 7*) w f 2dti k j - 2d 



E E 

3=1 aSA,|a 
2d-l 1 

/ 2d , E ^- 2d k^(2d-j)^ i e i/«i(« a /M a ) A 

3 : 

/2d,ifc" 



a(zA,|a|— j 



2d 

3=1 

2d— 1 

d E 

3 = 1 

< hd,i- 

Hence, (a Qji : a G A,i = 1, • • • , n) belongs to the feasible set of the geometric 
program in Corollarv l3.7l Plugging into the objective function, one sees after some 
effort that 



E(2rf-H) 



aeA 



2,1 



2d-l 

= E b i kj ^ ¥ ' 

3=1 



so r FK < f sos . 

Corollary 4.3. If \a\ < 2d for each a G A and fid.% > for i = 1, . . . , n, i/ien 

i 

^ i / , 2,! 

fgp — f dmt 



□ 



:=fo- E(2d-H) 



— I *' ' Q /2d 



2d 



where t :— |A|. 

Proof. Take a Qji = and apply Corollary 13.7 



□ 

Remark 4.4. Let C be a cone in a finite dimensional real vector space V. Let C° 
denote the interior of C. If a £ C° and b £ V then 6 G C° if and only if b— ea G C for 
some real e > (See [H Lemma 6.1.3] or 6] Remark 2.6]). Since £" =1 Xf d G E° dn 
O Corollary 2.5], for a polynomial / of degree 2d, with / 2c ; G E^d „, there exists an 
e > such that g = f 2 d ~ e(E"=i ^ 2<i ) e £ 2 d, n - The hypothesis of Corollary O 
holds for / — g. In this way, corollaries 13.71 |4~TI 14.21 and 14.31 provide lower bounds 
for f sos . Moreover, the lower bounds obtained in this way, using corollaries 14.11 ET21 
and 14.31 are exactly the lower bounds obtained in [6]. 



The bounds Tl, tfk, rdmt provided by corollaries 14.11 14721 and 14 . 3 1 are typically 
not as good as the bound f gp provided by Corollary 13.71 . 

Example 4.5. (Compare to [SJ Example 4.2]) 

(a) For f{X,Y) = X e + Y 6 + 7XY — 2X 2 + 7, we have r L « -1.124, r FK « -0.99, 
rdrnt w -1.67 and f sos = f gp » -0.4464, so f gp > r FK > r L > r dmt . 

(b) For f(X,Y) = X 6 + Y 6 + iXY + 10Y + 13, r L w -0.81, r FK w -0.93, 
r dmt ~ -0.69 and f gp « 0.15 « f sos , so f gp > r dmt > r L > r FK . 
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(c) For f(X,Y) = X i + F 4 + XY - X 2 - Y 2 + 1, / sos = f gp = r L = -0.125, 
tfk ~ -0.832 and r dmt « -0.875, so f gp = r L > r FK > r dmt . 
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